Protecting quantum correlations of two qubits in independent non-Markovian 

environments by bang-bang pulses 



o 

(N 

< 

(N 



-I— > ■ 



00 

m 

00 

o 



X 



Hang-Shi Xu and Jing-bo XrQ 

Zhejiang Institute of Modern Physics and Phystcs Department, 
Zhejiang University, Hangzhou 310027, People's Republic of China 

We investigate how to protect quantum correlations for two qubits each locally interacting with its 
own non-Markovian environment by making use of bang-bang pulses. It is shown that the quantum 
discord dynamics presents the phenomenon of sudden change for some certain initial states. We also 
find that the amount of quantum correlation between two qubits can be improved by applying a train 
of pulses and protected more effectively with shorter interval pulses or longer reservoir correlation 
time. 
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I. INTRODUCTION 

Quantum correlations has become an hot topic in 
quantum information and quantum computation and at- 
tracted much attention in recent years [1-5]. Entangle- 
ment is a special kind of quantum correlation and has 
been widely used in quantum information processing [1]. 
However, it has been reorganized that the quantum dis- 
cord which is firstly introduced by Ollivier and Zurck [2] , 
has an advantage over entanglement and can be respon- 
sible for the computational speedup for certain quantum 
tasks [3-6]. Quantum discord, which is defined as the dif- 
ference between the quantum mutual information and the 
classical correlation, is nonzero even for separate mixed 
states and may be regarded as a more general and funda- 
mental resource in the quantum information processing. 

Recently, the dynamics of entanglement and quantum 
discord for some open systems has attracted much at- 
tention [7-10]. In the non-Markovian environments, the 
memory stored in the environment has been considered, 
some of the initial quantum correlation which is lost dur- 
ing the dissipative dynamics can return to the qubits 
[7]. It has been shown that the discord is more ro- 
bust than entanglement under the Markovian environ- 
ments and vanishes only at some time points under non- 
Markovian environment [7, 8]. Furthermore, quantum 
discord has been predicted to show the sudden transi- 
tion between classical and quantum decoherence [9] and 
the phenomenon of sudden change [10] for some certain 
initial states. Some of these phenomena have been veri- 
fied by the recent experiments [10,11]. The enhancement 
of quantum discord for the system of cavity QED with- 
out dissipation by bang-bang pulses has been investigated 
[12]. Recently experimental suppressions of decoherence 
in a ring cavity and solid-state qubits using bang-bang 
decoupling technique have also been reported [13]. 

In this paper, we propose a scheme of protecting quan- 
tum correlations for two qubits each locally interacting 



with its own non-Markovian reservoir by means of bang- 
bang pulses. Using the master equation approach, we 
obtain an explicit expression for the reduced density ma- 
trix of the system. Firstly, we assume that the two qubits 
are initially prepared in the states with maximally mixed 
marginals [14] and show that the quantum discord of two 
qubits displays the phenomenon of sudden change for 
some certain initial states. Furthermore, we investigate 
how the bang-bang pulses affect the quantum discord and 
entanglement of two qubits and the phenomenon of sud- 
den change. It is found that the bang-bang pulses can be 
used to improve the quantum correlations of two qubits. 
The pulses with shorter interval can protect the quantum 
correlations of two qubits more effectively. The influence 
of reservoir correlation time and relaxation time scale are 
also investigated. 



II. DYNAMICS OF TWO QUBITS IN 
INDEPENDENT NON-MARKOVIAN 
ENVIRONMENTS WITH BANG-BANG PULSES 



In this section we investigate the dynamics of quantum 
correlations for two qubits each locally interacting with 
its own non-Markovian reservoirs with bang-bang pulses. 
The Hamiltonian of one qubit interacting with its own 
reservoir with bang-bang pulses is given by 



H = H + Hi + H P (t) 



(1) 



with 



Hn 



u k a\a k \ 



Ht 



E 

k 



where a k and a\ denote the annihilation and creation 
operators of the reservoir field and <j z — |e)(e| — |<?)(g|, 
cr + = |e)(g|, er_ = |g)(e| are the atomic operators. The 
index k labels different field modes of the reservoir with 
coupling constant g k and frequencies u> k . Hp is the 
Hamiltonian for a train of identical pulses of duration 



g k {a+a k + cr_4), 
(2) 
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r, i.e., 



H P {t) = Va z Y J 0{t-T-n{T + T))6{{n + l){T + T)-t), 

n=0 

(3) 

where T is the time interval between two consecutive 
pulses and the amplitude V of the control field is specified 
to be which means that we consider the 7r-pulse only. 
For the case that the pulses is strong enough, i.e. the 
duration r — > 0, the time evolution operator with pulse 
presented is given by Up — exp[— i(ff + Hj + t^(j z )t] ~ 
exp[-i|cr 2 ]. 

By making use of the master equation approach for 
the reduced density matrix, to the second order in per- 
turbation theory (Born approximation), we obtain the 
following equation for the reduced density matrix of the 
qubit [15] 



dp(t) 
dt 



= - [ <tfTr R [H I (t),[H I tf),ptf)®p R ]], (4) 
Jo 



where p(t), pp are the density matrix of the qubit and 
reservoir respectively, and Hj(t) is the Hamiltonian in 
the interaction picture which given by 



Hj(t) = U\t)HjU{t), 



(5) 



with U(t) = Te-'/o dt'[H +H P (t'))_ Usin g Eq . ( 3 ) ; we can 
obtain U(t) and Hj(t) as follows, 



U{t) = [U P U (T)] n Uo(t-nT), 



(6) 



H I {t)=Y,9k{-^P ] {e i ^ t a + a k + e- i ^ t a.al). (7) 



where n = [^] 

e -iH t 



denotes the integer part, and Uo(t) 



Assuming that the qubit is in resonance with the cav- 
ity mode, the reservoir spectral density is the Lorentzian 
spectral distribution [16]: 



J(u) 



1 



7oA 2 



2tt (wo - oj) 2 + A 2 ' 



(8) 



where 70 is the Markovian decay rate connected to the re- 
laxation time scale tr(~ Jq 1 ) and A is the spectral width 
of the coupling connected to the reservoir correlation 
time tb(— A -1 ). The Markovian and non-Markovian 
regimes are distinguished by 70 and A: 70 < A/ 2 means 
the Markovian regime and 70 > A/2 corresponds to the 
non- Markovian regime [16]. Using Eq. (4), (7) and (8), 
the master equation becomes 



dpjt) 
dt 



dt'K{t,t')Cp(t'), 



Cp{t) = (J-pa + - -cr+er-p 



1 



-pa+a- 



(9) 



(10) 



where C is the Markovian Liouvillian and K(t,t') is the 
memory kernel which given by, 



K(t,t') = (-l)^ ]+[t T ] k{t' -t), (11) 

where k(t) = j dujj^e 1 ^-^ 1 = ^e~ At is the corre- 
lation function. The factor (— is induced by the 
bang-bang pulses, which averaging to zero the unwanted 
interaction with the environment when T — > 0. 

The dynamics of a qubit S can be described by the 
reduced density matrix [15, 16] 



P 8 (t) 



(12) 



where the function P t obeys the differential equation 
dP /"* 

^ = j dt'K(t,t')P t ,. (13) 

The general solution of Pt can be derived as, 

xtr , d(t-nT) d(t-nT) l0 , N 
P t =e- xt lA n cos^^ l + Bn sm^^ >-} 2 , (14) 

where t <= [nT, (n + 1)T] and d = y/2j Q X - A 2 . Using the 
initial values A$ — 1 and Bo = ^ and boundary condi- 
tions, we can obtain recurrence relations of the constant 
coefficients A n and B n (n > 0), 

-Am , dT dT). 

A n = e 2 [A n _i cos— +S n _isin— J, 



2A 



e 2 [An-i cos — + sin 



dT 
T 



dT) 
2 J 



dT 



dT) 

^e"* J [A„_isin^- + B„_iCos-^-]. (15) 

In the standard product basis {|1) = |11),|2) = 
1 10), |3) = |01), |4) = |00)}, the elements for the reduced 
density matrix p AB {t) for the two-qubit system [16] can 
be constructed as, 

P?i B (t) = P? 1 B (0)P?, 
P22 B (t) = /4 S (0m + Pn(0)Pt(l P), 
p£f(t) = PifWPt + Pn(0)Pt(l Pt), 

pif(t) = 1 - [pffit) + P i 2 B {t) + p£ 3 B (t)}, 

Pi*(t) = P?2 B (0)P?,Pw(t) = P^(0)pf 
P?f(t) = P?4 B (0)PuP2 B (t) = (& B (P)Pt 

P24 B (t) = VPt[ P ^ B (0)Pt + P?f(0)(l Pt)}, 



Pt B {t) = VPt[pff(0)Pt + pf2 B (0)(l Pt)}, (17) 
and pf B (t)=pf B *(t). 



(16) 



t , 



3 



III. PROTECTING QUANTUM 
CORRELATIONS OF TWO QUBITS IN 
INDEPENDENT NON-MARKOVIAN 
ENVIRONMENTS BY BANG-BANG PULSES 

In this section we investigate how to protect quan- 
tum discord and entanglement for two qubits each locally 
interacting with its own non-Markovian reservoirs with 
bang-bang pulses. The definition of quantum discord is 
based on quantum mutual information which contains 
both classical and quantum correlations. For a bipartite 
system p AB , its total correlations can be measured by 
their quantum mutual information [2, 17] 



I(p A ")=S(p A ) + S(p*)-S(p A ») 



(18) 



where S(p) = —Tr{p\og 2 p) is the von Neumann entropy, 
and p A and p B denote the reduced density matrices of 
parts A and B, respectively. The quantum discord is 
defined as the difference between the quantum mutual 
information and the classical correlation [2, 17], 



Q(p A «) = I(p A «)-C(p A "), 



(19) 



where C(p AB ) is the classical correlation which depends 
on the maximal information obtained with measurement 
on one of the subsystems and can be expressed as 




FIG. 1: Qi(t) (solid line) and Q2(£)(dash line) are plotted as 
a function of time t with 70 = 1, A = 0.1 for different initial 
states (a)ci = 0.9, C2 = —0.9, C3 = 1 and (b)ci = 0.9, C2 = 
—0.9, C3 = 0.8 without control pulses. The analytical solution 
of the quantum discord is the minimum value assumed by the 
functions Qi(t) and Q2{t). 



d=l-P t + -{l + cz)Pl 



C(p AB ) = m^[S( P A ) - S{p AB \{B k }% (20) 

where {Bk} is a complete set of projectors preformed 
on subsystem B locally, S(p AB \{B k }) = Y.kPkS{Pk) 
is the quantum conditional entropy, pk = l/p k (I eg) 
B k )p AB (I ® B k ) is the conditional density operator and 
Pk = tr(AB)[(J <g> B k )p AB (I ® B k )] is the probability. 

We assume that the two qubits are initially in the 
states with maximally mixed marginals [14], which can 
be described by the three-parameter X-type density ma- 
trix 



P AB (0) = ~ 



1 + c 3 ci — c 2 

1 — C3 ci + c 2 

ci + c 2 1 — C3 

ci - c 2 1 + c 3 



(21) 



where c,(0 < c, < 1) are the real numbers. The state 
in Eq. (21) represents a considerable class of states in- 
cluding the Werner (|ci| = |c 2 | = | C3 1 = c) and Bell 
(|ci| = |c 2 | = I C3 ] = 1) basis states. Then, we obtain the 
density matrix p AB (t) for the two noninteracting qubits 
as 



p AB (t) = 



( a w 

b z 

z b 

w d 



(22) 



where 



-(l + C3 )P t 2 , 6=-P t --(l + c 3 )P t 2 , 



•z = 7(ci+c 2 )P t , w = 7 ( Cl - c 2 )P t . (23) 



Using the density matrix p AB (t), we obtain the explicit 
expression of quantum discord between two qubits as [7] 



Q(p AB ) =min{Qi,Q 2 }, 



(24) 



where 



h = S( P A ) - S(p AB ) - alog 2 (-^-) - Mog 2 (-^-) 

a + b a + 



"d + V 



d + b' 



(25) 



and 



Q 2 = S(p A ) - S(p AB ) - A+ log 2 A+ - A_; .g 2 A_, (26) 



-XL 



[cosf 
1\n-K 



with A± = \ ± ^(a-d) 2 +A(\z\ + \w\) 2 . 

If no control pulses is presented, Pt = 1 
■^sinY] 2 j which has discrete zeros at t r 
arctan(d/A)]/d, with n integer. It is easy to check that 
Q(t n ) — Qi(t n ) = Q2(t n ) = 0. For the time near these 
discrete zeros, for example at time t$ , we can obtain 
Qi(^) > Q 2 (£(7) by means of Taylor series expansion. 
For the initial state of Eq. (21), we have 

Qi(o) - Q 2 (o) = 1 + ( ~ 2 1)fc|c3 ' lo g 2 [i + (-i) fc |c 3 |] 

k=l 
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FIG. 2: Quantum mutual information 7(t)(a), classical cor- 
relation C(t)(b) and Quantum discord Q(t)(c) are plotted as 
a function of time t with 70 = 1, A = 0.1, ci = 0.9, C2 = 
-0.9, c 3 = 1 for different T: T = 0.4(dash line), T = 0.2 (dot 
line), T = 0.1 (dot-dash line) and without control pulses(solid 
line). 



-E 2 * log a [l + (-l)*x], (27) 
fc=i 

where x = max(|ci|, |ca|). We can check that Qi(0) < 
Qa(0) if M > max(|ci|,|c 2 |) and Q x (0) > Q a (0) if 
|c 3 | < max(|ci|, |c 2 |). If |c 3 | > max(|ci|, |c 2 |), noting 
that both and Qn(f) decrease monotonically when 

t < to as Fig- l( a ) shows, there is a point of intersec- 
tion of Qi(t) and Q^{t) when t is between and £0 ^ i-e., 
a peculiar dynamics of quantum correlation Q(i) with 
a sudden change occurs. However, no sudden change is 
observed if |c 3 | < max(|ci|, | ca j ) 9 as Fig. 1(b) shows. It 
should be noticed that, in this paper, the initial param- 
eters for the phenomenon of sudden change occur are 
different from Ref. [14] because of the different environ- 
ment. And this phenomenon is not found in Ref. [7] 
where the quantum discord dynamics of two qubits in 
independent non-Markovian reservoirs are investigated, 
but different initial states are chosen there. 

When the bang-bang control pulses field is presented, 
quantum mutual information, classical correlation and 
quantum discord are displayed as a function of the time 
t in Fig. 2. It is shown that the bang-bang pulses can 
protect both the quantum and classical correlations be- 
tween the two qubits. From Fig. 2(c), we can see that 
the time when sudden change occurs is prolonged by the 




FIG. 3: The concurrence C.E(i)(a) and the quantum discord 
Q(t)(b) of two atoms are plotted as a function of t and r for 
the parameters 70 = 1, A = 0.1 without control pulses. 
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FIG. 4: Quantum discord Q(t)(a) and classical correlation 
C(t)(b) are plotted as a function of time t with 70 = 1, A = 
0.1, r = 1 for different T: T = 0.6(dash line), T = 0.2 (dot 
line), T = 0.01 (dot-dash line) and without control pulses(solid 
line). 

bang-bang pulses. 

In order to quantify the entanglement dynamics of the 
qubits, we use Wootters concurrence Cs(t) [18] as a mea- 
sure. For the initial state of Eq. (21), the explicit expres- 
sion of the concurrence of two qubits can be obtained 
from the density matrix p AB (t), 

C E {t) = max{0, 2\w\ - 2\b\,2\z\ ~ 2Vad}. (28) 
Now, we choose c± — c 2 = c 3 = — r, i.e., the initial 
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FIG. 5: Quantum discord Q(t)(a) and concurrence Cf;(f)(b) 
are plotted as a function of time t with 70 = 1, A = 0.1, r = 0.5 
for different T: T = 0.6(dash line), T = 0.2(dot line), T = 
0.01(dot-dash line) and without control pulses(solid line). 

states of the two noninteracting qubits are the Werner 
state 

p w = r\ip-)(ilj-\ + ^I, (29) 

where \ip~) = (|0)a|1)b — \1)a\0)b)/V2 is a maximally 
entangled state and < r < 1. We plot the concur- 
rence and quantum discord as a function of t and r in 
the absence of control pulses field in Fig. 3. It can be 
seen from Fig. 3(a) that entanglement of two qubits is 
always zero for r < ^ [17] and vanishes and revives with 
time as r > |, which means that the entanglement sud- 
den death(ESD) phenomenon [16] appears for the sys- 
tem. However the quantum discord of two qubits van- 
ishes asymptotically as shown in Fig. 3(b). Different 
from Ref [12], here, both entanglement and quantum dis- 
cord cannot recover to the initial values according to the 
dissipation induced by multi-modes environment. How- 
ever, they can also be protected effectively by the bang- 
bang pulses with the same form. 

In Fig. 4. the quantum discord and concurrence as 
a function of t for different values of control pulses are 
displayed with the parameter r — 1, i.e., the two atoms 
are initially prepared in the maximally entangled state. 
We can see clearly that both the quantum and classical 
correlations between two qubits can be enhanced by the 
pulses with short time interval and the smaller interval 
between the pulses, the more enhancement of the quan- 
tum and classical correlations. For the Werner state with 
r = 0.5, the concurrence and quantum discord are plot- 
ted in Fig. 5 as a function of time t for different intervals 
of control pulses. We find that both the concurrence and 
quantum discord can be protected by applying the bang- 
bang control pulses. It is interesting to point out that the 
time when ESD occurs can be prolonged by the control 
pulses. 








FIG. 6: The quantum discord Q(t)(a) and classical correlation 
C(t)(b) are plotted as a function of time t and T with 70 = 
1, A = 0.1, r = 1. The functions is shown in the present of 
pulses(NeonColors surface) and in the absence of pulses(T- 
independent, LakeColors surface). 




FIG. 7: The difference of quantum discord AQ between Q(0) 
and Q(t) is plotted as a function of time t and (a) A~ x (~ tb) 
with 70 = 1,T = 0.01, r = 1 or (b) 7 ~ 1 (~ tr) with A = 
0.1, T — 0.01, r = 1. 
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We display how the quantum and classical correlations 
vary with time t and the time interval T in Fig. 6. It is 
shown that the pulses with short interval can protect the 
quantum and classical correlations between two qubits. 
Both the quantum and classical correlation descend al- 
most monotonously with time when bang-bang pulses are 
presented. When the interval T is long, the value of quan- 
tum or classical correlations would be less than the value 
in the absence of pulses for some region of time t. These 
phenomenon are similar to the quantum Zcno and anti- 
Zcno effects [19]. 

In order to explore the influence of the Markovian de- 
cay rate and the memory stored in the environment con- 
nected to the reservoir correlation time, we display how 
the difference of quantum discord AQ(t) = Q(0) — Q(t) 
changes as a function of A _1 (~ tb) or ■~f ~ 1 (— tr) in Fig. 
7. It is found AQ is smaller, i.e., the quantum correla- 
tions are protected more effectively by bang-bang pulses 
for longer reservoir correlation time tb or relaxation time 
scale tr. 

IV. CONCLUSIONS 

In this paper, we have proposed a scheme of protect- 
ing quantum correlations for two qubits in independent 
non-Markovian reservoir by making use of the bang-bang 
pulses. Considering the two qubits are initially in the 
states with maximally mixed marginals, we show that 



the quantum discord dynamics presents sudden changes 
for some given choices of initial parameters. We also ex- 
plore the influence of the bang-bang pulses on the quan- 
tum discord and entanglement of two qubits and the 
phenomenon of sudden change. The quantum correla- 
tions between two qubits are protected more effectively 
with shorter interval pulses or longer reservoir correla- 
tion time. Particularly, the time when ESD occurs can 
be prolonged by the control pulses and the time when 
sudden change occurs can also be prolonged. We only 
consider the ideal bang-bang pulses with the duration to 
be infinitely small and there is no interaction between 
two qubits. For coupled two qubits, the quantum corre- 
lations may be protected by the bang-bang pulses used 
in this paper. However, protecting the dynamics of cou- 
pled quantum systems is a problem needs further con- 
sideration since the decoupling pulses may disrupt the 
interqubit dynamics [20]. It is interesting to investigate 
the effect of finite duration time of the pulses and other 
pulse error for more realistic case [21]. The approach 
adopted here may be used to improve the implementa- 
tion of quantum information and quantum computation. 
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